Spin correlations in the paramagnetic phase and ring exchange in La2Cu04 
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Spin correlations in the paramagnetic phase of La2Cu04 have been studied using polarized neutron 
scattering, with two important results. First, the temperature dependence of the characteristic 
energy scale of the fluctuations and the amplitude of the neutron structure factor are shown to be 
in quantitative agreement with the predictions of the quantum non-linear sigma model. Secondly, 
comparison of a high-temperature series expansion of the equal-time spin correlations with the 
diffuse neutron intensity provides definitive experimental evidence for ring exchange. 
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Heisenberg was the first to realize that strong effec- 
tive spin interactions arise from the principle of the in- 
distinguishability of the particles Dirac generalized 
this concept, in the context of Group Theory, to include 
higher order interactions 2J. Multi-particle exchange 
dominates the physics of the quantum solid 3 HeQ but, 
surprisingly, is generally not taken into account for elec- 
tronic magnetic materials. The most powerful technique 
for exploring exchange interactions is the study of exci- 
tations from the ordered phase using inelastic neutron 
scattering. However, ambiguities in the interpretation of 
magnon dispersion curves sometimes mean that higher 
order terms remain hidden. By employing the indepen- 
dent approach of studying the instantaneous spin cor- 
relations in the paramagnetic phase, we obtain comple- 
mentary information that enables a better understanding 
of the exchange mechanism. In this letter we describe 
studies of the diffuse magnetic scattering from La 2 Cu04 
which provide compelling, quantitative evidence for the 
existence of four-particle cyclic exchange. 

La2Cu04 is of great intrinsic interest both as the 
parent compound of a canonical high-temperature su- 
perconductor, and as a very good realization of a 
two-dimensional quantum Heisenberg antifcrromagnet 
(2DQHAF). Magnetic Raman experiments^, infrared 
absorption studiesjp, and inelastic neutron scatter- 
ing measurements |3| show definitively the inadequacy of 
the nearest-neighbour Heisenberg model, and suggest the 
possibility that four-particle exchange may be significant. 

In an important series of experiments to study the 
diffuse magnetic scattering from La2Cu04 using un- 
polarized neutrons Q the temperature dependence of 
the magnetic correlation length was found to agree 
with the predictions of the quantum non-linear sigma 
model (QNLctM) 9J. However, the observed amplitude 
shows dramatic deviations from the predictions of this 
theory The QNLctM is the simplest possible effec- 
tive action for a 2DQHAF that is compatible with the 
long-wavelength spin waves and that does not assume 
a spontaneously broken symmetry. Moreover, its pre- 



dictions should hold even in the presence of four-spin 
exchange, as discussed below. Here we study the dynam- 
ical spin correlations above the Neel temperature using 
polarized neutrons, and find complete agreement of the 
observed 2D critical fluctuations with the predictions of 
the QNLctM. 

Dirac's approach provides the most transparent the- 
oretical framework to examine higher-order exchange 
interactions 2J. His analysis leads to an effective spin 
Hamiltonian 



Meff = -J2(- 1 ) PXJ ^° 
A 
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where A runs over all possible permutations of spin J"^ 
within the symmetric group, J\ is the exchange energy 
associated with a given permutation and p\ its parity. 
Any permutation can be expressed in terms of cyclic ex- 
change processes. Thouless was the first to point out 
that cyclic permutations of an even number of spins lead 
to AF exchange, whereas when an odd number of spins 
are permuted the resulting interaction is FM[lfJ. Only 
the values of the exchange parameters J\ depend on the 
choice of model; the form of interaction between spins is 
quite general. 

In La2Cu04, retaining the most important exchange 
processes involved in a plaquette, the general effective 
spin Hamiltonian is given by: 
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the J± > are pair-exchange frequencies between nearest 
(1), next-nearest (2) and next-next nearest neighbours 
(3), J3 and J4 represent three-and four-particle cyclic 
exchanges in a plaquette. In terms of spin operators 
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+(2J 2 (2) - 4J 3 + J 4 ) ^ Si • S, + 2J^ 3) ^ Si • Sj 
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We note that the four-particle cyclic exchange J4 in 
Eq. J5J contributes both four-spin and two-spin terms 
to Eq. J2Jl. On a square lattice, with the two-sublatticc 
antifcrromagnetic Neel phase there is a remarkable (al- 
though fortuitous) exact cancellation of all contributions 
of the J4 terms in linear spin-wave theory. This means 
that all quantities (in particular the magnon dispersion) 
calculated within this simple framework are the same as 
those corresponding to the pure Heisenberg Hamiltonian 
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with J^> = 4 X) -4J 3 , J£> = J^ J -2J 3 and J^°> = J y 2 
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and are completely blind to the four-particle permutation 
term J4. 

For simplicity we can model the Cu-0 planes in 
La2Cu04 using the half-filled one-band Hubbard model 
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where the hopping energy t characterizes the kinetic en- 
ergy, the potential energy U 3> t is the penalty for dou- 
ble occupancy, c (c + ) are the annihilation (creation) op- 
erators and n — c + c is a number operator. At fourth 
order in a k = t/U expansion, the J\s appear as 
4 1} /U = 2k 2 (1 + 4k 2 ), 4 2) /U = 12k 4 , 4 3) /U = 2k 4 , 
J 3 /U = 10k 4 and Ji/U = 20k 4 [HII3- More intricate 
expressions are obtained for a more general three-band 
Hubbard 

modelHU. 

The effective interaction between 
next-nearest neigbour pairs J 2 becomes negative (i.e. 
FM) because of the presence of the FM three-particle 
term J3. The next-next-nearest neighbour term is small 
and can be neglected. The magnitude of the four-particle 
cyclic exchange J4 is large. 

Since spin waves are insensitive to four-particle cyclic 
exchange, the curvature of the magnon dispersion at the 
zone boundary 7] is instead entirely due to the ferromag- 

~(2) 

nctic effective next- nearest neighbour exchange J 2 • hi 
contrast, there is no such cancellation of the four-spin 
term for the static susceptibility at high temperatures. 
We have, therefore, studied the diffuse scattering in the 
paramagnetic phase, and this is a new approach to the 
investigation of higher order exchange. The dynamical 
structure factor for neutron scattering is given by 



S(0) 



(cu-cq) 2 +T 2 (tu + cq) 2 +T 2 
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where T is the characteristic energy. Integration over 
energy transfer yields information on the equal-time spin- 
spin correlations since 

/oo 
S(Q,Lo)dcu » ^V^SfSg) « T X (Q). (7) 
-00 i 

The wave- vector dependent static susceptibility x(Q) carL 
be calculated from the exchange energies using a high- 
temperature series expansion. 

A 2g single-crystal of La2Cu04 from the array used to 
study the spin waves in the ordered phase was mounted 
inside furnaces, and the diffuse magnetic scattering was 
measured in the temperature range 300 — 500 K using 
the polarized neutron spectrometers D7 and IN20 at the 
Institut Laue-Langevin. XYZ polarization analysis was 
employed to separate the magnetic signal from the co- 
herent structural and spin-incoherent backgrounds|l4j . 
The scattering intensity measured in the (h, 0, 1) plane at 
room temperature using the multidetector on D7 is pre- 
sented in Fig.^a) for the nuclear scattering showing the 
structural Bragg reflections and (b) for the purely mag- 
netic signal showing the appearance of a rod of intensity 
perpendicular to the cuprate planes. The integrated in- 
tensity along the Q 3£) = (1,0, 1) rod for La2Cu04 was 
measured as a function of temperature with the incident 
wave vector fixed, ki — 2.08A" 1 , and the final wave vec- 
tor parallel to the normal to the cuprate planes in a sim- 
ilar manner to Ref.|| so that the cuprate square-lattice 
wave vector remained fixed at Q 2£) = (|, h ) for all energy 
transfers. For a quantitative temperature dependence of 
the intensity integrated over energy transfer it is essen- 
tial to determine how the spectral line shape varies with 
temperature. Energy scans were performed with Q fixed 
using the triple-axis spectrometer IN20 and typical spec- 
tra are presented in Fig. Eta). 

In the QNLctM, the correlation length £ is given by0 
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and the energy width T is related to £ by 
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where v s is the spin-wave velocity, p s is the spin stiff- 
ness and and Cr are undetermined constants of order 
one. Figure [5Jb) compares the temperature dependence 
of £ deduced from equation JJjJ using the values of T 
determined in energy scans. There is excellent agree- 
ment between the dynamical predictions of the QNLcrM 
and the correlation lengths measured using unpolarizcd 
neutrons by Birgeneau et al. 8J. The intensities mea- 
sured in fixed-Q energy scans on IN20 integrated over 



3 



energy transfer were converted to the amplitude S(0) 
in equation © using the known correlation lengths[8j. 
The intensities measured on D7 without energy anal- 
ysis were corrected using the spectral line shapes ex- 
trapolated from Fig.[2Ib), the instrumental energy win- 
dow, the Cu 2+ magnetic form factor and the correlation 
lengths, and the amplitudes from both experiments are 
combined in Fig. |2] The leading term in the expression 
for the ratio of the amplitude to the correlation length 
squared in the QNLuM is@ 
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FIG. 1: (Color online) Neutron scattering intensity in the 
(h,0,l) plane of La2Cu04 at room temperature, i.e. just 
above the Neel temperature, measured using the multi- 
detector on D7. Three-directional polarization analysis al- 
lows separation into (a) coherent structural scattering and 
(b) purely magnetic scattering, with removal of the incoher- 
ent background. A rod of magnetic scattering is developing 
along the [1,0, 1] direction showing the cross-over to 2D cor- 
relations. 



Figure l^a) shows that when data collected on D7 are 
corrected with a full knowledge of the spectral line shape 
they follow the same curve as those collected at fixed 
Q on IN20. Furthermore, the temperature dependence 
of intensities obtained using polarized neutrons is now 
in agreement with the predictions of the QNLcrM. The 
insensitivity of linear spin-wave theory to four-particle 
terms means that the expansion of four-spin exchange op- 
erators in terms of gradients of the Neel vector does not 



add any new terms to the QNLcrM. It is, therefore, grat- 
ifying that the clean measurements of the diffuse mag- 
netic signal using polarized neutrons agree now with the 
predictions of the QNLcrM for the renormalized-classical 
phase. 
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FIG. 2: (Color online) (a) Scans of energy transfer at fixed 
wave-vector transfer at several temperatures on IN20, and in 
(b) the characteristic energies are compared with the correla- 
tion lengths from Ref.0 using the QNLcrM^. 



In Fig. OUb) we show the comparison of the measured 
S(0) with the results of a high-temperature series ex- 
pansion of the multiple-spin exchange model Eq. J5J). 
The high-temperature series expansions were taken to 
fifth order and analytically continued using biased Pade 
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FIG. 3: (Color online) The neutron scattering amplitude 5(0) 
at Q = (|, |) of the cuprate square lattice, (a) The ratio of 
5(0) over the square of the magnetic correlation length || 
varies linearly with temperature squared, in agreement with 
the QNLo-MQ. (b) The temperature dependence of 5(0) fol- 
lows a straight (solid) line and the gradient agrees with the 
static susceptibility calculated in high-temperature series ex- 
pansions with J a, — 0.25 Jj 1 ' (dashed line). The dash-dot line 
shows the calculation with J4 = 0. Fourth and fifth order 
expansions are shown, the latter extending to lower tempera- 
ture. 



approximants|lq| . The values of the pair exchange ener- 
gies are those corresponding to the effective pair exchange 



2J. 



(i) 



111.8 meV and 2J 2 (2) = 



-11.4 mcV deduced from 



the magnon spectrum in the ordered phasejg. Plotted 
on a semi log scale to extract the leading behaviour in 
1/T, the experimental results fall on a straight line (solid 
line), and the gradient is in perfect agreement with the 
predictions of the series expansion at high temperatures 
with J4 = 0.25J 2 1 ' ) (dashed line) derived from the Hub- 
bard model. This agreement is achieved with no free 
parameters except an overall scale factor. The dramatic 
difference in slope with respect to the dash-dot theoreti- 
cal line (obtained with J4 = 0) demonstrates the extreme 



sensitivity of the diffuse magnetic scattering to this term. 
These data constitute the first quantitative evidence for 
four-spin cyclic exchange in La 2 Cu04. We note that the 
ratio J4/J2 « 0.25 is compatible with the predictions 
of the one-band Hubbard model, but more accurate neu- 
tron data would allow comparison with a more general 
three-band model|l3j|. 

The higher order terms found to be of crucial im- 
portance in the physics of solid 3 He are also shown to 
be significant in an electronic magnetic material. It 
seems highly likely that ring exchange will be impor- 
tant in many other electronic magnetic systems, espe- 
cially in those with strong hybridisation paths, such 
as the CU4O4 plaquettes. Optical experiments indicate 
that higher order exchange is important in other high- 
temperature superconductors including YBa 2 Cu30fi 
Bi 2 Sr 2 Cao. 5 Yo. 5 Cu 2 8+y , Nd 2 Cu0 4 and Pr 2 Cu0 4 [ 
The magnitude of the four-spin cyclic exchange is com- 
parable to the pairing energies, and it is possible that cir- 
culating electronic currents have an important role in the 
mechanism of superconductivity. Ring exchange is also 
believed to be important in related ladder compounds, 
such as La6CagCu 24 04i and SriACu^Ovi [lq. 

In summary, four-spin cyclic exchange has been re- 
solved in diffuse scattering experiments in the paramag- 
netic phase of La 2 Cu04, and the 2D critical fluctuations 
are correctly described by the QNLcrM. 
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